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Abstract 



It is well-known that the LIE(Locally Induction Equation) admit soliton-type solutions and 
same soliton solutions arise from different and apparently irrelevant physical models. By com- 
paring the solitons of LIE and Killing magnetic geodesies, we observe that these solitons are 
essentially decided by two families of isometries of the domain and the target space respectively. 
With this insight, we propose the new concept of geometric solitons of Hamiltonian flows on 
manifolds, such as geometric Schrodinger flows and KdV flows for maps. Moreover, we give 
several examples of geometric solitons of the Schrodinger flow and geometric KdV flow, includ- 
ing magnetic curves as geometric Schrodinger solitons and explicit geometric KdV solitons on 
surfaces of revolution. 
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1 Introduction 



The dynamics of vortex filaments has provided for almost a century one of the most interesting 
connections between differential geometry and soliton equations, and an example in which knot- 
ted curves arise as solutions of differential equations possessing an infinite family of symmetries 
and remarkably rich geometrical structures. The first mathematical model of the evolution of an 
Eulerian vortex filament, an approximately one-dimensional region of the fluid where the velocity 
distribution has a rotational component, was derived in 1906 by Luigi Da Rios [6], a student of 
Tullio Levi-Civita. Both from mathematical and physical points of view, one has found that the 
locally induction equation (also called Da Rios equation) is a fascinating nonlinear evolution equa- 
tion. It is related to the well-known Landau-Lifshitz (LL) equation which describes the dynamics 
of several important physical systems such as ferromagnets, moving space curves, etc. and has 
intimate connections with many of the well known integrable soliton equations, including nonlinear 
Schrodinger and modified KdV equations. It admits rich dynamical structures including solitons, 
dromions, vortices, etc. 

A vortex filament is a curve 7(i, s) in the Euclidean space M 3 depending on time t and arc-length 
parament s£R. The evolution of 7 in Locally Induced Approximation (LIA) is determined by the 
following Da Rios equation which is also called Locally Induction Equation (LIE): 

dj 57 <9 2 7 , „ 

m=^ x d^> or 74 = 7 X7 ' 

A careful derivation of the self-induced dynamics of a vortex filament can be found in G. K. 
Batchelor [2j. H. Hasimoto [2] discovered in 1972 that the LIE is equivalent to a complex- valued 
cubic Schrodinger equation via the so-called Hasimoto transformation: 



$(7) = kexp J rds 



where k is the curvature and r is the torsion of 7. 

It is well-known that the cubic Schrodinger equation defined on R admits a traveling wave 
soliton of the form 

$(i,s) = $o(s-c£). (1.1) 

From the correspondence between the solutions of cubic Schrodinger equation and the solutions 
of LIE, it is not difficult to see that (jl.ip gives a soliton type solution 7(3, t) to the LIE whose 
curvature k and torsion r satisfy 

k(t,s) = k(s - ct), r(t,s) = t(s - ct). (1.2) 

Inserting (jl.2p back into the cubic Schrodinger equation via Hasimoto transformation, we reduce 
the LIE to an intrinsic soliton equation of k and r. 

However, this intrinsic equation is not quite illuminating and the geometric significance is not 
obvious. Moreover, it is always difficult to recover the curve 7 from the intrinsic equation of k 
and r. Hence, a natural problem is how to understand the soliton equation from the viewpoint of 
global differential geometry and reveal the hidden geometric structure of such solitons. We observe 
that the connection between solitons of LIE and Killing magnetic geodesies sheds some light on 
this problem. 
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The so-called magnetic geodesic is just the trajectory of a charged particle moving on a Rie- 
mannian manifold with a magnetic field. The electron experiences a Lorentz force F( ■ ) and the 
magnetic geodesic 7 satisfies equation 

Vy 7 ' = F(Y). 

In particular, if the manifold is three dimensional, then the Lorentz force is related to a vector field 
V and the equation of magnetic geodesic is 

V y7 / = FA 7 / , (1.3) 

where A is the cross product induced by the Riemannian metric. If the vector V is a Killing field, 
then a solution to the above equation is called a Killing magnetic geodesic. 

Many mathematicians and physicists have been interested in the existence of closed geodesies 
on a manifolds. It seems that to study the existence of closed magnetic geodesies on a manifold is 
a more challenging problem. Even for the case of two dimensional manifolds, it is a hard and long- 
standing problem and related to the problem of finding curves with prescribed geodesic curvature. 
There is a large amount of works on this field using different approaches in the last decades, and the 
results on existence of closed magnetic geodesies turn out to be subtle under various assumptions. 
We refer to [211 \TT\ [25] and references therein. 

On the other hand, the LIE can also be well-defined on three dimensional manifolds by 

7 i = 7' A Vy Y. (1.4) 

In [3], the authors investigated the Killing magnetic geodesic in three dimensional space forms 
carefully and by verifying the intrinsic equations, they found that the Killing magnetic geodesies 
are just solitons of the LIE. However, this coincidence seems somewhat surprising and mysterious. 
In fact, one has also found that some other important but apparently irrelevant physical models, 
such as elastic curves, share the same reduced equation for their respective solitons. We refer to 
|10| for more details and examples. 

In this paper, our primary goal is to provide a more global and geometric view on the connection 
between solitons of LIE and magnetic geodesies. The key observation is that the soliton of LIE is 
essentially decided by two families of isometries. Namely, they are of the following form: 

j(s,t) = (fi t 070 oip t (s), (1.5) 

where (fit is a one-parameter group of isometries of the ambient space, Vt( s ) * s a one-parameter 
group of isometries of the domain, and 70 is a curve independent of time t. 

If we assume that (fit is generated by a Killing vector field V, then by putting the form (II. 5p 
back into the LIE (jl.4p . we will get a reduced equation for 70, which is exactly equation (jl.3p . In 
the specific setting of classical LIA, the ambient space is the Euclidean space M 3 and (fit consists of 
translations and rotations, while the domain is K and 

ip t (s) = s - ct 

is just a shifting of the parameter s with speed c. Indeed, the soliton 7 is generated by rotating 
and translating the initial curve 70, which turns out to be a Killing magnetic curve, with a slipping 
motion along the curve itself. 
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On the other hand, the LIE can be viewed as a Hamiltonian flow and there is an important gen- 
eralization, namely, the Schrodinger flows for maps from a Riemannian manifold into a symplectic 
manifold O [T5l [30] . The Schrodinger flow is an important geometric flow which is closely related 
to the harmonic map and wave map. However, there seems to be no analog of soliton solution to 
the Schrodinger flow except, in some literature, the term "soliton" refers to harmonic maps which is 
just a stationary solution of the flow equation. The significance of our observation is that a solution 
of the form f)1.5|) essentially gives a unified definition of solitons for geometric flows on manifolds. 

More precisely, suppose (M, g) is a Riemannian manifold and (N, uj) is a symplectic manifold 
with a compatible almost complex structure J. The Schrodinger flow for a map u : [0,T)xM->JV 
is 

— = J(u)t(u), 

where 

r(u) := tracegVdu 

is the tension field. Then the soliton type solution of the Schrodinger flow analogous to (jl.5p should 
has the following form: 

u(t,x) = (j) t °u oipt(x), (1.6) 
where tpt and 4>t are one-parameter groups of isometries of M and ./V respectively, and 

is a map independent of time t. Since this kind of soliton solution is closely involved with the ge- 
ometry of the underlying manifolds, we should call them "geometric solitons" . (See Definition 14.11 
below for more details.) We will see that many important and special models arising in mathemat- 
ical physics fall into the scale of geometric solitons, which provide us solid evidences to justify the 
name and significance of geometric soliton. In particular, we found that the geometric soliton of 
1-dimensional Schrodinger flow on Riemannian surfaces turns out to be a magnetic geodesic, and 
the well-known equivariant solutions of the Schrodinger flow for maps from R 2 or S 2 into a Kahler 
manifold with axial symmetry (see [12] ) has the same structure as geometric solitons. 

The geometric soliton can also be defined for the geometric KdV(Korteweg-de Vries) flow which 
is recently proposed by Sun and Wang [29J. The geometric KdV flow also origins from the movement 
of vortex filaments, this time along with an axial flow, which transforms to a classical mKdV equa- 
tion under the Hasimoto transformation. From the viewpoint of integrable systems, the geometric 
KdV flow from 1R into § 2 belongs to the same hierarchy of the corresponding Schrodinger flow. 
Similarly, the geometric KdV soliton is defined as a solution of the geometric KdV flow in the form 
(|1.6p . In fact, this kind of geometric solitons can be defined for all similar Hamiltonian flows on 
manifolds. 

The next goal is of course seeking the above geometric solitons. By inserting (II. 6p back into 
the flow equation, we are led to study an equation of the generating map uq which is independent 
of the time variable t. For the Schrodinger flow, we get an elliptic system and uq turns out to be 
a harmonic map with potential. For the geometric KdV flow, we get a third-order o.d.e. equation. 
Each solution no to the reduced equation generates a geometric soliton of the geometric flow. 
However, in general the reduced equation is very hard to solve except some special cases where the 
underlying manifolds are simple. 

This paper is organized as follows. First we recall the classical traveling wave solutions of the 
LIE of vortex filament with or without axial flow, as well as elastic curves as solitons of LIE. Next 
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we show the equivalence of these solitons and Killing magnetic geodesies extrinsically using our 
observation fjl.51) . Then after a brief review of the of Schrodinger flow and geometric KdV flow, 
we move on to give the definition of geometric solitons. Meanwhile, we derive a reduced equation 
for the geometric solitons. Finally, we show some examples when the underlying manifolds are 
simple. This includes magnetic geodesies as geometric Schrodinger solitons, and explicit examples 
of geometric KdV solitons on surfaces of revolution. 

2 Solitons of LIE 

2.1 Localized Induction Approximation 

The motion of a thin vortex filament in an inviscid incompressible fluid is described asymptotically 
by the so-called Locally Induction Approximation(LIA). This is a minimum model which was first 
discovered by Da Rios in 1906: 

It = i X V'. (2-1) 

If 7 is parameterized by arclength and k, r denotes the curvature and torsion of 7 respectively, 
then under the Frenet frame {t,n, b}, we have 

7 ' = t, t' = kn, n = -kt + rb. 

Thus equation (|2.ip becomes 

It = kb. 

Remark 2.1. If s is the arclength parameter 0/7(0, •) at the initial time t = 0, it is easy to verify 
that s is arclength parameter of "y(t, •) for all time t > 0. 

Hasimoto found a famous transformation which is now referred as the Hasimoto transformation 

$(7) = fcexp (i J rds'^j . (2.2) 

Using this transformation, equation (|2.ip collapses to a cubic Schrodinger equation for the complex 
valued function <I>: 

-m = ^ ss + ^\ 2 + A)<i>, (2.3) 

where A is a real function of time variable t (which can by eliminated by shifting the origin). 

Note that the function $ only depends on the intrinsic quantities of the curve 7. Namely, <J> is 
given by the curvature k and torsion t which is independent of the choice of coordinates and the 
parameter of 7. Thus equation (12. 3p should be regarded as the intrinsic version of the LIE (|2.ip . 

A remarkable feature of the cubic Schrodinger equation (|2.3|) is that it admits a type of solitary 
wave solutions which demonstrate a kind of self-similarity. These solutions have the form 

$(t,s) = $o(s-ci) (2.4) 

where c is a constant describing the speed of the wave and 3>o satisfies an elliptic equation: 

ic$(, = ^ + i(|$ | 2 + ^o. 
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Obviously, such solutions correspond to a type of solitary wave solution j(t, s) of the LIE (|2.1|) with 

k(t,s) = k(s - ct), r(t,s) = t(s - ct). (2.5) 
Putting this back to equation (|2.3p , we get the intrinsic soliton equation for k and t as the following: 



k" + h 3 - kr 2 + A(t)k = -ck(r - r(-ct)), 
2k' t + kr' = ck'. 

In the next subsection, we will see that the intrinsic equation of elastic curves coincides with the 
above equation. 

2.2 Elastic curves as solitons of LIE 

An interesting and remarkable discovery noted by Hasimoto in [14] is that the soliton of LIE is 
related to elastic curves. The problem of elastic curve was first proposed by Bernoulli to Euler in 
1744. In mathematics, the bending energy of a thin inextensible wire 7 in 3 dimensional Euclidean 
space is described by the total squared curvature 



E(j) = J k 2 ds, 



where k is the curvature of the curve 7. A curve is called an elastic curve if it minimizes the 
bending energy E. 

In practice, we often use a modified energy which depends on the differentials of the curve 
explicitly. Namely, we define an elastic curve to be a critical point of the following energy: 

F{l) = \ I W'\ 2 ds+ I K{s)(W\ 2 -l)ds, 



2 _ 

where A(s) is a point-wise Lagrangian multiplier, constraining the speed of the curve. Obviously, 
such an elastic curve in equilibrium satisfies the Euler-Lagrangian equation: 

7 ""-l(A7') = 0. 
ds 

Suppose 7 is parameterized by arc-length and using the Frenet frame along the curve, the above 
equation can by rewrite as: 

k" + - A; 3 - kr 2 - -k = 0, 

2 2 (2.7) 

2k' t + kr' = 0. 

Here A is a constant which is determined by the Lagrangian multiplier 

A = — k 2 + -. 
2 2 

The second equation integrates to give a point-wise constant quantity 

k 2 T = c, 
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and the first equation becomes 

fcSfc" + 1 fc 6 _ ^4 _ c 2 = (2 g) 

With the help of elliptic functions, one can solve equation (|2.8|) to get (refer to [21)] ) 

(!-£■«(£.,,)), (2.9) 

where sn(-,p) denotes the elliptic sine function with parameter p, and the parameters p, w and &o 
are related to constants A and c satisfying < p < w < 1 and 

J 2Aw 2 = A;g(3u; 2 - p 2 - 1), 
l4eV = ^(l-™ 2 )(™ 2 -/). 



Now let us recall the intrinsic equation (12.6P for solitons of the LIE. If we introduce a new 
variable £ := s — ct which is independent of i. Then the only two terms dependent of t in the first 
equation of (|2.6ft should cancel with each other. Namely, we have 

A(t)k - ckr(-ct) = A x k, 

where A\ is a constant. Comparing equation (|2.6[) to equation (|2.7p . one can easily find that a 
soliton of the LIE with c = is exactly an elastic curve. Thus we get a special family of solitons of 
the LIE given by ®M . 



Remark 2.2. Langer and Singer give a complete classification of elastic curves in M 3 . The elastic 
curve may also be defined on 3 dimensional manifolds and similar results can be obtained. Refer to 
126)1 and references therein. Thus we may obtain many potential solitons of the LIE. However, it 
is somehow difficult to recover all the formations of possible solitons from the expression of k and 
t. This was accomplished by Kida in fit 



2.3 Vortex filament with axial flow 

A vortex tube is often accompanied with an axial flow in the core. This phenomenon was modeled 
by Moore and Saffman [20] and then simplified by Fukumoto and Miyazaki [9j as the following 
equation: 

It = «7' x 7 " + PW" + \l" x (V x t")L (2-10) 

where a and (3 are two constant numbers. In fact, the right hand of the above equation is a linear 
combination of the first two vector fields of the Localized Induction Hierarchy [17J. Under the 
Frenet frame, equation (|2.10p reads: 

7t = akb + p(^k 2 t + k'n + kTb), (2.11) 

It's obvious that when a = 1, f3 = 0, equation (|2.10p reduces to the LIE (|2.ip . On the other 
hand, if a = 0, /3 = 1 the equation becomes: 

7i = 7 "' + ^ 7 "x( 7 'x 7 "), (2.12) 
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Applying the Hasimoto transformation (|2.2p . equation (|2.12p can be rewrite as a complex valued 
mKdV equation: 

$t = $ sss + ~|$| 2 $ s . (2.13) 

This is the second equation in the nonlinear Schrodinger hierarchy and it also admits a solitary 
wave soliton which satisfies (|2.4p . and hence (|2.5p . Following the same spirit, one can compute the 
intrinsic equation of a soliton type solution to equation ()2. 12[) : 

3k" t + 3k' V + kr" - kr 3 + -k 3 r = -ck(r - r(-ct)), 

3 , , (2 ' 14) 

fc'" — 3k! t 2 — 3/ctt' + ijk 2 k' = — cfc'. 



3 Killing Magnetic geodesies 
3.1 Killing magnetic geodesic 

A magnetic field on a Riemannian manifold (M,g) is characterized by a closed 2-form $7 € A 2 (M) 
which denotes the set of all smooth section of A 2 (T*M). A moving particle of charge in the presence 
of a magnetic field F will experience a Lorentz force F which is represented by a skew-symmetric 
(1, 1) tensor defined by 

g(F(X), Y) = n(X, Y), VX, Y e T(TM). 

Here r(TM) denotes the set of smooth vector fields on M. 

The magnetic geodesic 7 describes the motion of an electron in a magnetic field which satisfies 
equation 

VyV = F(7 / ). (3.1) 

In particular, if the manifold M is 3 dimensional, then the 2-form VI is equivalent to a vector 
field V via the Riemannian metric g. Indeed, one can define V = (*f2)*, where * is the Hodge star 
operator and * denotes the promotion given by the metric g. Moreover, dVl = implies 

div(V) = trace 9 (W) = 0. 

On the other hand, we can also define a cross product using the volume form by the formula: 

g (x ay,z) = n (x,Y,z), yx,Y,z G r(TM). 

In this case, it is easy to check that the Lorentz force is F(X) = V A X. Thus the equation of 
magnetic geodesic (|3.ip becomes: 

VyV = V A y. (3.2) 

Especially if the vector field is a Killing field, that is, V is a vector field which generates a family 
of 1-parameter group of isometries, then a solution 7 to equation (13. ip is called a Killing magnetic 
geodesic. Note that the Killing magnetic curves can be explicitly determined, see [18 ^ [23" 1 [2^]. 

In [1], the authors discovered a Killing magnetic geodesic with arc- length parameter in 3 space 
form satisfies the following equivalent equation of f|3 . 2|> : 

k" + -k 3 - kr 2 + Ak = -ukr, 

2 (3.3) 

2k' t + kr' = uk', 
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where A is a constant number and u is the so-called quasislope defined by 



u) = (V, t) = const. 

It is easy to see that this intrinsic equation (|3.3|) is exactly the same with the equation derived 
for the solitary wave solution (|2.6f) of LIE. The interesting connection between vortex filament and 
Killing magnetic geodesic in three-dimensional space forms has been addressed in [3]. But their 
observation is base on comparing the intrinsic equations of the curvature and the torsion, which 
seems somewhat obscure and mysterious. In the next subsection we are intended to explain the 
equivalence of Killing magnetic curves and solitons of LIE from a more global and extrinsic view 
point, which seems more natural and illuminating. This will motivate the definition of geometric 
soliton in the next section. 

3.2 Relation to solitons of LIE 

Now we try to explain the equivalence of Killing magnetic curves in Euclidean space M 3 and solitons 
of LIE. Suppose 7(s, t) is a solitary wave solution of the LIE which comes from the soliton of the 
corresponding Schrbdinger equation via the Hasimoto transformation. Then the curvature and 
torsion of 7 satisfies equation (|2.5p . i.e. 



k(s,t) = k(s — ct), r(s,t) = t(s — ct). 

We need to derive the underlying equation for the soliton 7. Note that the above equation does 
not necessarily imply that the curve also satisfies 

j(s,t) = j(s - ct). 

Indeed, both the curvature k and the torsion r are intrinsic properties of the curve. Namely, they 
stay unchanged if the curve 7 varies by an isometry, which corresponds to a translation and a 
rotation in the Euclidean space M 3 . Therefore, we may suppose that there exists a one-parameter 
group of isometries 4>t and a curve 70 such that 

7(5, t) = M-Yoi* ~ <*))■ (3.4) 

If the one-parameter group of isometries <pt is generated by a Killing vector field V, then we may 
compute 

j t l(s,t) = V(Mlo)) + #t(7o)(-c7(>) (3 5) 



On the other hand, we have 



i (F( 7o )) + ^(7o)(-c7o)- 



^7M) = #*(7o)K) (3.6) 



and 



'2 



d 

-r^7(s,t) = dMloKlo) + d 2 <M7o)(7o,7o) 



ds 2 

Since 4>t is an isometries, we have d^cptilo)^, 7o) = 0. Hence 



^7(^,t) = # i (7o)(7o')- (3-7) 
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Putting ([53]) . (|3T5T> and ([32]) back to the LIE flST), we can get 



d«M7o) x dMO = dMV(7o)) + #t(-c 7 o) 



which implies 



7o x To = V(tq) - c 7o . 



(3.8) 



Comparing equation (|3.8p with equation (|3.2p . we find out that they are equivalent. That is, if 7 
is a solitary wave solution of the LIE f)2 . If) . 70 turns out to be exactly a Killing magnetic geodesic 
in M 3 . Or inversely, every Killing magnetic geodesic 70 in R 3 generates a soliton type solution 7 to 
the LIE. 

4 Geometric Solitons 

4.1 Geometric Hamiltonian flows 

It is well-known that there exists a direct relation between the Landau-Lifshitz equation and the 
Da Rios equation. One can obtain the Landau-Lifshitz equation by differentiating the Da Rios 
equation if the vortex filament curve is parameterized by arc-length parameter. About a century 
after the discovery of the Da Rios equation, mathematicians began to realize that the Landau- 
Lifshitz equations stemming from vortex filaments can be generalized to a Hamilton flow for maps 
from a Riemannian manifold into a symplectic manifold, i.e., the Schrodinger map flows for maps 
into a symplectic manifold. 

Recently, Sun and Wang [29] proposed another geometric flow now referred as the geometric 
KdV flow which belongs to the same hierarchy of the Schrodinger flow when the domain manifold 
is M and target is a compact Hermitian symmetric space [30]. For the reader's convenience, we 
include here the definitions of these flows. 

Suppose (M, g) is a Riemannian manifold and (N, oj) is a Symplectic manifold. Assume N 
is embedded in an Euclidean space W K , J is a tamed almost complex structure and h is the 
corresponding Riemannian metric such that 



Then the tangent space of W at a map u € W is the pull-back tangent bundle T U W = u*(TN). 
There is a natural symplectic form Q on W induced by u). Namely, for any X, Y G T U W, 



as a Hamiltonian function on W, then the associated Hamiltonian vector field Xe 1 satisfies 



u>(t,r)) = h(Z,Jr)),VZ,T)€T v N. 



Let W be the Sobolev space defined by 



W := {u G W 1 ' 2 (M,'R K )\u{x) G N for a.e. x G M}. 




Consider the energy functional 




Q{X El ,-) = dE 1 (-). 
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Thus we can define the Hamiltonian flow by 

u t = X El (u) = JVE 1 (u). (4.1) 

It is well-known that 

Tl (u) ■= VEi{u) = Au + A(u)(X7u, Vu) 

is the tension field where A(u) denotes the second fundamental form. So the Hamiltonian flow (14.ip 
of Ei is just the Schrodinger flow 

u t = Jn(u). (4.2) 

When the dimension of M is 1, namely, when M = R 1 or M = S , there is another important 
functional, namely, the pseudo-helicity, given by: 

E 2{u) = - (\7 x u x ,Ju x ). 

1 J M 

Similarly, one can also define the Hamiltonian flow of E<i by 

u t = Jt 2 (u), (4.3) 

where 

r 2 (u) := VE 2 {u) = -JV 2 x u x + -R(u x , Ju x )u x 

and R is the curvature tensor of N. This is the geometric KdV flow defined by Wang and Sun 
recently, see [29J for more details. 

Next let's explain the relation between the above geometric flows and LIE. Suppose "f(x) is a 
closed curve in M 3 and is parameterized by arc-length, then u := 7' defines a map form S 1 into S 2 
since |n| = = 1. Recall that there is the Frenet frame {t,n, b} along 7 in IR 3 with t = 7'. On 
the other hand, there is an orthonormal basis {T, N} on S 2 along u where T is the unit tangent 
vector along u and N is the normal vector. Note that the complex structure on § 2 is J := ux and 
we may suppose N := JT. Then, we have 

u x = 7" = t = kn = vT, 

where k is the curvature of 7 and v := \u x \ is the speed of u. Thus 

n = T, k = v. (4.4) 

Furthermore, 

V x u x = V(„t)(i>T) = v x T + v 2 k g N 

= ( 7 '") T = (k'n - k 2 t + krb) T 
= k'n + krh, 

where k g is the geodesic curvature of u. Without loss of generality, we may assume that k ^ 0, 
then 

b = N, k g = ~. (4.5) 
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Now one can easily find that in this case, the Hamiltonian functions are just integrals of the 
curve 7 that arise from the study of elastic curves. Indeed, 



E\{u) = \ \u x \ 2 ds = / k 2 ds. 
Js 1 Js 1 

E2(u) = / (X7 x u x ,u x u x )ds = I k 2 rds. 
Js 1 Js 1 

On the other hand, after differentiating with respect to x, the Da Rias equation (|2.ip becomes 

Ut = U X u xx 

which is just the Schrodinger flow (|4.2p on S 2 . Moreover, the LIE of vortex filament with axial flow 
(|2,l(jp becomes 

3 

u t = aux u xx + (3[u xxx + -(u x x (ux u x )) x ]. 

One may easily check that the latter part of the above equation is exactly the geometric KdV flow 
on S 2 . 

Remark 4.1. The geometric flow naturally connects the two apparently irrelevant but important 
physics models: elastic curves and vortex filaments. This relation remains valid in the case of 
three dimensional space forms, see Jfflj . The remarkable coincidence was stressed by many authors, 
but usually observed by comparing the intrinsic equations. Our perspective via the geometric flows 
provides a more global and geometric view. Namely, the LIA hierarchy is the Hamiltonian flow 
corresponds the bending energies of elastic curves and rods. 



4.2 Geometric solitons 

Inspired by the relation between solitons of vortex filaments and Killing magnetic geodesies, we 
propose a new concept of geometric soliton of the Hamiltonian flows, which can be regarded as a 
geometric generalization of the classical solitary wave solitons. 

Definition 4.1. Suppose (M,g) is a Riemannian manifold and (N,h,J) is a Kdhler manifold. A 
solution u to the Hamiltonian flow Q or j4-3\ ) is called a geometric soliton, if there exist two 
one-parameter groups of isometrics ipt : M — >■ M , (fit '■ N — > N and a map v : M —> N independent 
of t such that u has the form: 

u(t,x) = (j) t O V o ip t (x). 

One can also define a geometric soliton which only involves the isometries of the target manifold. 
Namely, solitons which have the form u(t,x) = (fit v(x) where (fit is a family of isometries of the 
target manifold. This kind of geometric solitons of the Schrodinger flow was first noted by Ding 
and Yin [8] in order to find periodic solutions to the Schrodinger flow into S 2 . Song and Wang 
|28j investigated geometric solitons of Schrodinger flow into Lorentzian manifolds and successfully 
obtained soliton solutions of the hyperbolic Ishimori system. It is worthy to point out that the 
existence of geometric solitons is conjectured to be related to the first eigenvalue of the underlying 
manifolds. 
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The concept of geometric solitons is a natural extension of the solitary wave solution to the 
LIE found by Hasimoto, or the classical solitons of the Schrodinger equation and KdV equation. 
Indeed, a solitary wave solution of the LIE satisfies (|3.4p . i.e. 

70,*) = <M7o(s - ct)), 

where <pt is the isometry of the 3 dimensional Euclidean space consisting of translations and rota- 
tions. The only difference at first glance is that the traveling speed c of the solitary wave is missing 
in our definition of geometric solitons. The effect of c corresponds to a slipping motion along the 
curve itself. However, if we regard the transformation 

ipt{s) = s - ct 

as an isometries of the real line, the consistence of the definition of geometric soliton and classical 
solitary waves is justified. On the other hand, for a linear Schrodinger equation 

iut = Au 

defined on a flat torus T m , a solitary wave solution is of the form u(t) = e lkt v where A; is a positive 
constant and v is a real function which satisfies the equation 

Av + kv = 0. 

Here, <f>t = e tkt can be viewed as an isometric group and again this classical solitary wave solution 
falls in the scale of above defined geometric soliton. 

Another important example is the well-known equivariant solution of Schrodinger flow on man- 
ifolds with axial symmetry. The m-equivariant map from M 2 to S 2 has the form 

u(r,9) = e mm v(r), 

where (r, 0) are polar coordinates on M?, m is an integer, v : [0, oo) — > S 2 is a map indepen- 
dent of t, and R is the matrix generating rotations around the u 3 -axis. In a given homotopy 
class, the m-equivariant harmonic map Q m minimizes the Dirichlet energy. A general problem on 
Schrodinger flow which has attracted a considerable attention for the past ten years is to describe 
the flow for initial data near Q m , see for example [12J. A major break-through on Schrodinger flow 
recently made in [19J is the discovery of a finite time blow up profile. The formation of singularity 
corresponds to the concentration of a universal bubble given by 

u (x,t) = e^ R Q 1 ^y (4.6) 

If we regard <pt = e Q ® R and ipt(x) = jfm a s one-parameter group of transformations, then the 
bubble (|4.6p has the same structure with our geometric soliton. 

After modulo the isometries of the underlying manifolds, the only thing left is the time- 
independent map v. Just as the relation of solitons and Killing magnetic curves shown in the 
last section, the map v satisfies a reduced equation. This is formulated in the following theorem. 
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Theorem 4.1. Suppose the one-parameter groups of isometries {(fit}, {^t} ore generated by a holo- 
morphic Killing vector field V € T(TN) and a Killing vector field W £ T(TM) respectively. Then 

i) u = (fit ° v o tfi t is a geometric soliton of the Schro dinger flow {-j.Sfy iff v satisfies 

Jn(v) = V(v) + V w v, (4.7) 

ii) u = (fit o v o ifit is a geometric soliton of the KdV flow f^.3| ) iff v satisfies 

Jt 2 {v) = V{v) + V w v. (4.8) 

Proof. Since the proof of the above two equalities are essentially the same, here we only prove i) 
for the geometric Schrodinger soliton and ii) follows. 
The geometric Schrodinger soliton u has the form 

u(t,x) = (fit ov o ifitix), 

where (fit and ifit are one-parameter groups of isometries generated by V and W respectively, that 
is, 

^( y ) = v((fi t (y)), Vyeiv 

and 

^-(x) = W{Sh{x)), Vi£M. 
at 

The group property of (fit implies 

(fit O (fis = (fis+t = <fit° <fis- 

Differentiating the above equality with respect to s at s = gives 

d(fi t {V) = V((fi t ). 

Similarly, 

dMW) = W{ifi t ). 
Hence by differentiating u with respect to we have 

du 

— (t, x) = V{(fi t ovo ifi t ( x )) + d(fi t o dv{W{ip t {x))) 
ot 

= d(fi t (V(v(rfit{x)))) + d(fi t (V w v(Mx)))- 

On the other hand, for any t and x := ifit(x), we have 

V 2 u(t,x) = V(d(fi t (Vv(x))) 

= d(fi t (V 2 v(x)) + d 2 (fi t {Vv(x), Vv(x)). 

Since (fit is an isometries, the second fundamental form of (fit is zero and it is obvious that the latter 
part of the above equality vanishes. Thus by taking trace, we get 

Tt(u(t,x)) = d(fit{T\{y{x))). 

Recall the assumption that the family of isometries (fit is holomorphic, which means Jod(fit = d(fit°J- 
Consequently, the Schrodinger flow ut = Jt\(u) now becomes 

d(fi t {V{v{x)) + V w v(x)) = J o d(fi t (n(v(x))) = d(fi t (Jn(v(x))), 

which is equivalent to 

Jt\(v) = V(v) + Vh/d. 

□ 
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5 Examples 



5.1 Geometric Schrodinger solitons on surfaces 

Here we discuss the special case of one dimensional geometric Schrodinger solitons on surfaces. 
Suppose (M 2 ,g) is a Riemann surface and 

u = <p t o 7 o ip t : R 1 x [0, T) -> M 2 

is a geometric Schrodinger soliton, where <pt and tpt are one-parameter groups of isometries and 7 is 
a curve with arc-length parameter s. If (f>t is generated by a Killing vector field V and ipt = s — ct, 
then by Theorem 14.11 the generating curve 7 satisfies the reduced equation: 

J(7)V s7 , = V( 7 ) - c 7 „ (5.1) 

where J is the complex structure of M 2 . Since 7 is parameterized by arc- length, we may denote 

T:= 7s , N:=JT and V s j s = k g N, 

where k g is the geodesic curvature of 7. Then from the equation (|5.ip we may deduce that there 
exists a function g such that V(7) = 9(7) 7s- Therefore equation ()5.ip becomes 

V s7s = (<?(7) - c) J( 7 )7»- (5-2) 

The above equation shows that the curve 7 is a magnetic geodesic on M 2 . In fact, on a 
Riemannian surface M 2 , a magnetic field Q. has the form Q, = kQo where Oo is the volume form of 
M 2 and k is a function. Thus the corresponding Lorentz force F is F(X) = kJX and the equation 
(|3.ip for magnetic geodesic becomes 

V s7 , = fc(7)A7)7«- (5-3) 

Comparing equation (|5.2j) with (I5.3p . we may find the Schrodinger soliton 7 is just a magnetic 
geodesic with 

Hi) = 9(7) - c - 

It is obvious that a magnetic curve on M 2 with arc-length parameter has geodesic curvature 
k. Thus magnetic geodesic is closely related to the problem of prescribed geodesic curvature. In 
particular, a geometric Schrodinger soliton with trivial Killing vector field V = is a curve with 
constant geodesic curvature c. 

It has been a long-standing problem in the last decades to find closed magnetic geodesies 
or curves with prescribed geodesic curvature on surfaces which was initiated by Arnold [TJ. The 
existence turns out to be subtle and somehow related to the topology and geometry of the underlying 
manifolds. Various results are obtained by using different approaches, including Morse-Novikov 
theory and symplectic methods, etc. For example, Schnieder [25] successfully applied an infinite- 
dimensional index theorem to prove that there exist at least two simple closed magnetic curves 
on S 2 if certain geometric assumptions are assumed. This would certainly yield the existence of 
two geometric Schrodinger solitons on S 2 . For more results on closed magnetic curves, we refer to 
[111 [2~T| [25] and references therein. 
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5.2 Geometric KdV solitons on surfaces of revolution 

Next we deduce the equation of geometric KdV soliton surfaces of revolution in the three-dimensional 
Euclidean space R 3 . 

To begin with, let £ be a surface of revolution in R 3 with coordinates z±,Z2 and £3. Without 
loss of generality, assume that £ is obtained by revolving a regular plane curve 7(7*) in the plane 
Z2 = around the Z3-axis with r as its arc-length parameter, say (f(r),0,g(r)). Then £ could be 
expressed as 

£ =: (/(r)cos0,/(r)sin0,0(r)), 

with 

/ r 2 + 5r 2 = 1, 

where f r , g r are the derivatives of f(r) and g(r) respect to r respectively. 

Let <f>t(z) : 1 x S -> S be the obvious one-parameter group of isometries on £ which is just a 
rotation around the Z3-axis with speed oj given by 



M*) = 

Then the corresponding Killing vector field on £ is 

V(p) = u)f(— sin 0, cos 6, 0) T , p = (/(r) cos 0, /(r) sin 9, g(r)) € £. 

Moreover, let ipt(x) = x + ct be the one-parameter groups of isometries of R with constant speed c, 
then the corresponding Killing field W is just the constant c. By Definition 14, 1\ a geometric KdV 
soliton u(x, t) : R x [0, +00) — > £ on £ has the form u(x, t) = o u o tp t where u = u(x) : R — > £ 
is the initial map. For simplicity, we consider the 1-equivariant map u with 9 = x and r = r{x). 
Then by Theorem 14.11 the initial map 

u{x) = (f(r(x))cosx,f(r(x))smx,g(r(x))), 

satisfies the reduce equation 

^l u x + 2 R ( u xi J u x)Ju x = V(u) + V w u- 

Here Vwu = cu x and 

u x = (r' f r cos x — f sin x, r'f r sin x + f cos x,r'g r ), 

where r' , 9' denotes the derivative respect to x. Since £ is a surface, it is easy to check that 

R(u x ,Ju x )Ju x = G(u)\u x \ 2 u x , 

where G{u) is the Gauss curvature of u(x) on £. Note that on the surface of revolution £ the 
Gauss curvature 

G(u) = G(f) = 
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only depends on /, which means G is constant along a parallel. Now we may focus on the equation 

G{u) 



vx + 



\u x \ u x = V(u) + cu a 



(5.4) 



Let v := \u x \ be the speed of u(x), T = t^t be the unit tangent vector of u(x) and N € TE the 
unit tangent vector orthogonal to T. Then we have 



v 



u, 



r' 2 + f, 



and 



V x u x = V x {vT) = v'T + v 2 k g N; 



V 2 x u x = V x (v'T + v 2 k g N) 



= (v" - v 3 k 2 )T + (3vv'k g + v 2 k' g )N, 
where k g is the geodesic curvature. Here we used the formula 

V T T = k g N and V T N = -k g T. 
Substituting (f53|) into yields 



{v" - v 3 k 2 + 



G(f), 



Let 



c)T + (3vvk g + v 2 k' g )N 



V(u) = k\T + k 2 N, 




it is easy to check that 



h = V(u) ■ T 



k 2 = V(u) ■ N 



iofr' 



V V 

In view of (|5.7p and (|5,6p we get the equations of geometric KdV soliton as follows 



// 4 7 2 , "a/) 4 f 2 

vv — v k„ H — v — cv = uf ; 

y 2 

3v 2 v'k g + v 3 k' g = -ufr'; 
v = r + j . 



We are interested in the following special cases: 

i) If the KdV soliton has constant speed, i.e. v = vq, then (|5.8|) becomes 



4 7 2 i "C/J 4 i-2 

' — — v - cv = uf ■ 



V Q K g 



r' 2 + f 2 



(5.5) 



(5.6) 



(5.7) 



(5.8) 
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ii) If (ftt = id, i.e. the Killing field V = and u = 0, then the geometric KdV soliton is just a 
KdV curve studied in [27] and 

' vv" - v A k 2 g + -cv = 0; 

< y 2 

k 3v 2 v 'k g + v 3 k' g = 0. 
Thus (v s kg)' = which implies i> 3 fc 9 = C is a constant. Thus we have 

Moreover, in this case, if u(x) has a constant speed, i.e., v = vq is constant, then the geodesic 
curvature on the KdV curve must be constant and so is Gauss curvature. In another word, the 
curve should be a parallel of £ with r(x) = r . 

hi) Inversely if we assume the KdV soliton is a parallel of S, i.e. r(x) = r$ is a constant, then 
both / = /(ro) and G = G(f) are constants. The equation (j5.8|) yields 

// 4 7 2,^4 c2 

i;?; — v k + —v — cv = uj ; 
y 2 

* 3vVA; 9 + v 3 k' g = 0; 

2 j?2 
V = f . 

which implies that v = f is constant. By the same arguments, the quantity t> 3 /c s = C is constant 
and so is the geodesic curvature k g . Therefore equation (|5.8p reduces to an equation of /(ro) and 
f"(r ): 

\f"f+u;f + cf + C 2 = 0. 

Thus given a surface of revolution and a constant C, we can use the above equation to identify a 
parallel which is a geometric KdV soliton with v 3 k g = C. 

5.3 Geometric KdV solitons on cylinders 

In order to see the geometric KdV soliton more precisely, we now show explicit examples on a 
cylinder. Suppose r > is a constant and 

S = (r cos 9, r sin 6,h), 6 6 [0, 2tt], h G R 

is a cylinder of radius r. There are two natural Killing fields on the cylinder generated by the 
rotation and the translation along the Z3-axis respectively. Namely, a Killing field on S is the 
combination of a rotation and a translation and has the form: 

V = uR + aT 

where R = (— sin 6, cos 9, 0) T , T = (0, 0, 1) T and u>, a € R. The corresponding one-parameter group 
of isometries (j>t(z) : R x S — > S is given by 

(cos cut — sino;t \ / z\ \ I z\ \ 
sinwt coswt I £2 I + ^2 
1 J \ z 3 J \z 3 + at J 
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Let u{x) = (r cos kx, r sin kx, h{x)), x G R be a curve on £, where k 6 M is some constant. The 
metric on a cylinder is flat and the curvature vanishes. Thus the equation of geometric KdV soliton 
on S is just 

Uxxx = V(u), 



or more precisely, 



This implies uj = rk 3 and h" 



rk 3 sin x 

—rk 3 cos x 
•in 





f — sin x ^ 






H 


cos X 


M 






\ o I 







h'"{x) 

x) = g. Hence we obtain that 



a 



h(x) = '-x 3 + C 3 x 2 + C 2 x + Ci, 


where Ci{i = 1, 2, 3) could be any constants. 

If a = 0, then we obtain rotating solitons. In this case, h(x) is a linear combination of x 2 , x and 
1. Obviously, if h(x) = C%, the curves are just parallels on the cylinder; if h(x) = C 2 x, (C2 7^ 0), 
then we get helixes; while if h(x) = C3X 2 , (C 2 7^ 0), the geometric solitons are parabolic curves 
rolled on the cylinder. If a ^ 0, the one-parameter isometric groups contains a translation along 
the 213-axis. These KdV solitons are cubic curves on the cylinder. In fact, the rotating KdV solitons 
are combinations of these three kinds of curves. (See Figure 1) 






(a) h(x) = C 2 x (b) h(x) = C 3 x 2 (c) h(x) = C 4 x 3 

Figure 1: Geometric KdV Solitons on Cylinder. 
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